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Abstract 

■  The  g-tenSors  (2"**  order)  of  doublet  radical  gallium  arsenide  clusters  Ga^ASy  (x  +  y  =  3,  5),  were 
calculated  from  first  principles  using  new  packages.  Geometry  optimizations  and  hyperfine coupling  constants 
are  also  r^orted,  using  the  B3LYP/6-3IlH-g(2df)  level  ESR  results  were  compared  to  btperimental  Ga^ASj 
data,  and  previous  calculations  for  GajAs,  Ga  As2,-Ga2As3.  New  ESR  and  structural'results  are  presented  for 
GaAs4  and  Ga4As.  Our  results- for  GajAsj,  the  only  Ga-As  cluster-  far  wiiich  experimental  ESR  data  exists  were 
a  better  fit  to  experimental  values  thanprevioas  calculations,  implying  our  Ag  results  for  the  other  Ga^ASy 
clusters  are  also  of  good  quality. 
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1.  Introduction.  .  . 

To  date  there  have  been  fewejtperimental  ESR  studies  on  ni-V  doublet  radicals,  namely  only  on  BNB' 
and  GajAsj  ^  to  our  knowledge.  For  triplet  or  quartet  ni-V  radicals,  experimental  ESR  data  have  only  been 
fround  fcM-  Gai As*  (X^S'),’  GaP*  (X^S),*  GaAs- AI,;Ga,.,As,’  anrisite  drfects  (To, in  GaP,‘  Asg,  in  GaAs,’  and 
P,„  iniiP*),  andGa-vacancies  in  electrcn-irradatedGaP.’  • ’  •  ;  •  ‘  '  : 

However  there  have  been  a  number  of  theoretical  studies  on  the  electronic  states  and  structures  of  Ill-rV  . 
radicals,  in  particular  for  Ga-As  containing  neutral  and  ionic  doublet  radicals. Arratia-P6rez  and 
Hemandez-Acevedo'^  calculated  the  magnetic  Zeeman  and  hyperfine  interactions  of  GajAsj  using  the  self- 
consistent  Dirac  scattered  wave  method  (SCF-DSW-Xa)  of  Yang  et  aF  and  a  first-order  perturbation 
•  procedure;  which  coiifirm^  Weltiier’s^  ESR  spectral  determination  of  a  trigonal  bipyramidal  sthicture  for" 
GaijAsj;  Arratia-Perez  and  Hera^dez-Acevedo  havft  his  6  calculated  the  inagnetic  Zeeman  and  hyperfine' 
interactions  for  GaAsj  and  Ga2As,'^althou^  to  datethefe  areno  experimental  ESR  data  foF these  species. '  • 

We  compare  these  calculated  and  expa  imental  data  to  our  results,  and  also  use  the  g-teiisor  results 
•tas  a  launching  point  to  the  study  of  larger,  clusters  which  approach  bulk  properties,  with  and  without  site 
defects.  .  •;  '  ■  '  '.  '  '  •'  '  '  ■  • 

2.  Methods  ’  . 

Geometry  optimizations  and  hyperfine  coupling  constant  (hfee)  calculations  were  carried  out  with  tlie 
GAUSSIAN  98  suite  of  programs^*  at  the  B3L^'P/S-3  1 1+G(2df)  level.  Starting  geometries  were  those  niven 
in  the  literature  for  Ga  As,,'''  '’  Ga, As, Ga..^s,  and  GajAs,,'’-’®  GaAsj  and  Ga^As,'’  with  other  possible 
geometries  also  e.xamined.  For  GaAsj  and  Gaj.A.s.  only  skeiehes  orthe.struciures.wcre  provided.'' 

The  ihcorelical  evalualion  ofg-lensors  using  periiiibation  tivory  is  described  in  detail  in  rel's.  26  and 
.27.  Th.e  .total  Ag  fora  given  molecule  is  comprised  of  first-  and  secondrordertehns.  fn  this, paper,  only  seepnd- 
order  g-lensor  components  were  ealeuhited,  as  the  first-order  eonlributions  to  the  total  Ag  tire  e.xitecled  to  be 
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small  (ca.  lOQ  ppm)  in  rdation  Id  the  second  older  Ag  values.  The  contribution  to  Ag  (2"'')  is.diie  to  the 
coupling  of  an  excited  state  with  the  ground  state,  and  is  proportional  to  their  spin-orbit  coupling  (SO)  and 
magnetic  transition  moment  (L)  matrix  elements,  and  inversely  proportional  to  their  energy  separation  (AE). 
The  total  second-order  Ag.is  calculated'as  a  sum-ovef-stafes  expansion,  which  generally  involves  strong 
•  coupling  to  only  the  first  few  low-^hg  excited  states.**  ’" 

.  Programs  developed  earlier.in  this  group*’  worked  very  well  for  P'  and  2"*  order  g-tensor  calculations 
(e.g.  refs .  27  -  30)  but  were  limited  to  systems  with  up  to  ca.  80  electrons.  Ga-As  containing  clusters,  the  focus 
of  this  paper,  quickly  surpassed  this  limitation.  Also,  speed  is  an  issue  if  the  calculation  of  larger  clusters  is 
to  be  practical.  New  programs  were  acquired  and  modified  for  our  purposes  to  calculate  the  g-tensors  of  these 
Gai-As  clusters.'  T^ey  are  based  oti:  the  Turbomole  package’’  for  efficient  integral  and  SCF  calculations;  qn 
the  GrimmemultifeferenceCI(MRCI)package,’*  which  also  gives  the  a^ular  momentum  (L), matrix  elefnentS' ' 
that  we  require;  aiid  finally  on  the  Marian-Hess  mean-field' spin-orbit  integrals”  as  implemented  •  by 
.  Schimmelpfennig,”  and  adapted  for  tlie  Grkmne  MRCIpacka^  by  Marian  and  Kleiaschmidt.”  Here  the  one- 
and  two-electron  spin-orbit  elements  are  caiculated  from  an  effective  one-electron. one-center  mean-field 
approximation.  Therefore  the  spin-orbit  matrix  elements,  which  in  fuH  cfh  initio  mode  require  large  amounts 
of  computer  time,  can  be  calculated  in  much  less  time  but  are  not  as  accurate.  .  ■  ■  ,  ■ 

The  electronic  charge  centroid  (EC C)  was  taken  as  gauge  origin.”’  '  ' 

,  •  Test  results  with  the  new  package  gave  g-shifts  within  15-20  %  of  MRCI”  values  for  small  molecules 
containing  fir.st-  and  second-row  atom.s  (e.g.  ’A10,  H,C0C  LiCA.”'  With  the  new  package,  the  following  TTT-V 
doublet  radicals  were  .studied:  GaAs.  (97‘electroas),  Ga.As  (95  e).  Ga^As,  (Ibl’e),  Ga,As.  (159  e),'GaA.Sj 
( 1 63  e),  and  GujAs  ( 1 57  o'). .  ,  ..... 


3.  Opliniized  geumeti  k's 


••  •  •  ■  .  .  ■  .  .  .  ...  ,  ..  •  •  .  5 

.  All  geometry  optimfeations  were  dohQ  at  the  B3LYP/6-3 1 1  +g(2df)  level  of  theoiy.  The  results  for  the 
lowest  energy  structures  are  given  in  Table  1,  and  compared  with  those  of  previous  calculations.  The  various 
structures  are  shown  in  Figure  1. 

•  3.1.  GaAs2  •  ■  •  •  ■  .  . .  .  . 

GaAsj  was  first  examined  in  1987,^®-and  again  in  2000.”  by  Balasubramanian  as  a  triangular  Q, 
structure;  a  1991  stucfy  also  ccnsida-ed  a  linear  geometry.'*  Work  by  Meier  et  al  in  1991'“'  examined 
triangular  C,,  and  linear  C„,. )  geometries.  In  all  cases,  the  ground  state  was  ^B,  with  a  Cj,  triangular 
geometiy. 

Our  results  also  gave  in  Cj,  symmetry  as  thd  lowest  energy  state.  Alternate  possible  geometries 
considered  here  for  GaAsywere  linear  (X’llj;  As-Ga-As;  Ga-As  =  2'.22  ‘A)  ahd  C„v  (X^IIVGa-As-As:  Ga- ' 
As  =  2.48.A,  As-As  =  2. 16  A)  structvfres  that  were  1 .62  and  0.52  eV,  respectively,  higher  in  energy  tlan 
(Figure  1).  This  canpares  well  to  tlie  results  of  Meier  e/  o/,''*  where  X’Bj  vyas  lower  than  tire  (2.39  A)  and 
C„„  (Ga^As  =  2.61  A,  As-As '=  2.24  A)  structures  by  1.53  and  0.81  eV;  respectively.- 

3.2.  Ga^As  '  •  '  .  •  .  •  ■  ' 

For  Ga^.As.  Balasubramanian  reported  an  X^B,  ground  state  with  C,,'  symmetry  and  an  angle  of  79.9°,' 
witha  C|.  (A!)  structure  (Ga-As  =  2.283,  2-534  A;Ga-As-Ga  =  ■90.3°)'  0.025  eV  Irigher  in  ener^  (MRCI  + 
Q,  whe-eOwas  defined  as  multireference  Davidson  cOirections  tothe'MRCT  energies  for  uncoupledquadaiple  ' 
clusters).”  Two  other  low'lying  states,  ^B,  (2.52  A,  108.2°,  AF,  =  +0.22  eV)  and  -  A,  (2.47  A  1 18.5°,  AE  = 
+0.I9.OV.).  wcrealso  rqiortai  inref  17.  Bahisubramanian  Ibunda  ‘B,  ground  .stale  in  his  1991  study."' 

.  •  ..  Our  .altenpts  to  .r.qno.duce  Balagubramanian’s  C,,  (X'B,)  geometry  horn  ref  1:7  .(2.0Q0)  resulted,,, 
instoki  in  an  (Ja-.-Xs-Cia angle  o\' cci.  96  tlegrees  |  I'or  the  6-3 1 1  •  g(2dr)  basis  set,  B3  LX'P  gave  Ga-.'Xs  2.384  A, 
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96.3°;  MPW1PW91  gave2.358  A,  93.3°;  and  MP2  gave  2.365  A,  103.5°]  for  GajAs,  inclose  agreement  with 
the  local  spin  density  (LSD)  result  of  98.0°  from  ref.  19  (see  Table  1).  Our  linear  D.^  (X^n„,  2.447  A)  and 
C„^,  (X'll,  As‘Ga  =  2.246  A,  Ga-Ga  =  2.665  A)  structures  were  0.09  and  0.84  eV  higher  in  energy, 
respectively,  than  our  X^Bj  Cjv'bent' geometry.  InBalasubramaniah’s  1991  paper  with  me  *B,  ground  state,’* 

■  ..the  ordering  was  Q,  (X^B,:  2.527.A,  109.5°)  0  eV;  D-^  (2.498  A)+0.12  eV;  C,,,  (As-Ga  -  2.319  A,  Ga-Ga 
•  =2.735  A)-K)..68eV.  ...  .  .  .  .....  ... 

iJ.  Ga.Asji 

The  structure  of  Ga2As3  was  precictedin  1992  to  be  a  trigonal  bipyramid  with  an  X^A^"  ground 
state,  using  the  LSD  method.*^  The  ESR  spectrum  \vas  obtained  one  year  later,  from  which’ a  trigonal' 
bip  jo'amidal  structtirewas  profKised.  ^  MRCI  calculations  by  Liao  et  alf^  arid  HF  followed  by  MP  2  calculations  . 
by  Piquini  et:al,^^  also  resulted  in  trigonal  bipyramidal  structures.  * 

Our  calculations  showed  the  lowest  energy  structure  of  Ga2As3  to  have  an  X“ A^"  ground  state  with 
trigonal  bipyramidal  geometry,  in  agreement  with  the  st ructu  re  propos  ed  from  the  ©cperimental  ESR  data^  and  . 
previous  calculations.*^'’^’^^- Alternate  possible  gcometries-consideredby  us,  shown  in  Figure-!,  Were:  a  C,, 
-square pyramid  97  eV  higher  Ilian  the  structure;  and  a  Q  square  pyramid  (5^ A”,  two.G^  in  the  , 

‘base’)  0.46  eV'higher  in  energy.-  .  ’  •*'  ’  ‘  ‘ 

3f4/Gn,As. 

CASSCF  and  MRGT  calculatiais  by  using  relati\istic  effective  core  potentials  (RECP), 

sliowed  two  nearly  degenerate  isomers  for  Ga,As2:  a  distorted  trigonal  .bipyramid  C.,  (a)  in  Figure 

1 ),  and  a  C,,  edge  capped  tetrahedron  (XMl,,  C.,  (b)  in  Figure  I),  with  the  trigonal  bipyramid  being  lower 
in  ena’gy  by  0.005  eV  (MRCI  j-  Q,  and. 0.03  eV  CASSCF).  LSD  calculations  by  Lou e/ a/ gave  similar  results,  . . 
with  an  energy  sejxiration  of  0.01  eV.*'* 
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For  GajAsj,  we  also  found  that  the  lowest  energy  ^metry  was  not  as  for  Ga^As,,  but  rather  a 
Cjv  dislortedtrigonal  b  pyramid  havii^  a  ^A,  ground  state  with  the  equatorial  Ga  atoms  arranged  in  an  isosceles 
triangle;  the  energy  gap  between  the  and  Dj^  geomefries  was  0. 1 3  eV  (0. 1  eV  in  ref.  19).  Other  geometries 
we  investigated  (see  Figure  1)  werea  Cj,  edge^capped  fetrahedron  [^B,,  Cj,  (b)  in  Figure  1]  6 16  eV  higher 
in  energy;  a  Q  square  pyramid  <X?A”,  two  As  in  the  ‘base’)  0.64  eV  higher  in  energy;  and  a  Cj^  square  ' 
•  pyrsi'uid  [X^A,,  Cj,  (c)  in  Figure  l-J,  1 .36-  eV  higher  than  the.  Cjy  distorted  trigonal  bipyramid  [Cjy  (a)]. 

5.5.  GaAs^ 

GaAsj  was  included  in  a  study  of  electronic  and  structural  trends  in  small  GaAs  clusters  by  Piquini 
et  dl,  who  carried  but  HartreeTpocIc  S.CF  calculations  including  all  electrons  and  tio  symmetry  constrairits 
duritig  bptirhization,  followed  by'  siri^e-point  MP2  calbulatiofis  bn  the  minimurn  energy  ccrifigurations.'’ 
However,  they  did  not  report  any'geomelry  detail's'other  than  a  sketch  of  the  molecule,  the  symmetry,  and  a 
table  of  selected  bond  orders. 

•  Our  results  gave  a  Q;  edgchcapped  tetrahedron  Gj;  (a)  in  Figure  1]  as  the  lowest  energy 

structure,  resembling  that  shown  in  rbf.  15.  We  calculated  square  pyramidal  (J^A,)  and  Cj/plartar 

trapezoidal  [X.'B,,  Cj,  (b)  in  Figure  1]  geometries  to  be0.30  and  1 .68  eV,  respectiydy,  higher  inenergy  tlian 
the  Cj,,  edge-capped  tetrahedron.  '  ‘  '  •  •  •  ’  ' 

5.6.  GcijAs  '  '  '■ 

GujAs  was  also  included  in  the  stuck'  by  Piquini  et  a1'^  and  again  no  detailed  geometrical  information  ' 
was  givai.  I  he  lowest  energy-structure  we  obtained  for  Ga,.-\s  was  a  C;>  planar  trapezoidal  sti'ucture  {X' )  . 
with  the  As  atom  in  the  center  |C,,  (a)  in  Figure  1  ],  similar  to  the  structure  obtained  by  Piquini  ei  cil.'-  Our 

results  showed  square  pyramidal  .(Q„  X.-B|)-and  C2,.edge-.capped  .tetrahedral- [X^AzGCj.  (b)-in  Figure  1.] 
geometries  to  be  respectively  0.23  and  0.56  e\'  higher  in  energy  than  the  planar  irapezoidti  1  struciure. 
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4.  ^ypetfine  muplii^  constants . 

The  atomic  charges,  spin  densities  (SD)  and  hfcc  data  are  given  in  Table2.  For  all  six  molecules,  the 
Ga  atoms  carry  positive  charges  and  all  As  atoms  are*  negatively  charged,  as  expected  since  As  is  more 
electronegative  than  Ga.  The  SOMOsare’shown  ih  Figure  2.  .  ‘  * 

4.L  GaAs2  ,  ■  . .  -  .  .  .  . .  '  ,  •  .  ,  ...  . 

The  SOMO  of  GaAs.  is  composed  of  an  in-phase  combination  of  a  Gap^  orbital  with  As  p^  orbitals, 
fonning  a-type  antibonding  orbitals  along  each  Ga-As  bond,  as  also  reported  in  ref.  13.  These  Ga-As 
antibonding  orbitals  of  the  SOMO  areca.  parallel  to  the  r-axis,  perhaps  due  to  tlie  small  As-As  distance  (2. 193 
A)  which  is  slightly  shorter  than'that  calciilated’for  As;^  (2'.30  A)  knd  As2'  (2.31  A).‘^  ‘  ‘ 

'  The  SD  is  grfeater  on  die  As  atoms  (0:368  per  atom)  than  on  the  Ga  atom  (0.262);  This  is  ih  gbbd 
agreement  with  the  spin  populations  calculated  in  ref.  13.,  who  report  0.338  (0. 106  isotropic,  0.232  anisotropic) 
per  Ga  and  0.331  (0. 121  isotropic,  0.210  anisotropic)  per  As  atom.. The  anisotropic  spin-dipolar  couplings 
(Txx.yy.zz  Tablc  2)  are  much,  larger  than  Ai^^  for  both  Ga  and  As,  and  the  relative  magnitudes  of  our  spin-  ■ 
dipolar  couplings  are.  similar  to  those  . from  ref:  13  .(see -Table  2), 

4.2.  GaAs  .  '  '  ■  ■  ’  • 

*  ■  The  SOMO  is  formed  from  an  in-phase  combination  of  a  p,  orbital  from  all  three  atoms  according  to 
ref.  \3,  but  this  corresponds  to  the  SOMO-l  in  our  calculations:  oiif  SOMO  is  comprised  of  three  in-plane  p- 
orbitals:  p^  on*As  and  the  p,,  orbitals  bfGa  lie  along  the  A*s-Ga  bonds  such  that  tlie  negative  lobes  for  all  three 
orbitals  point  toward  the  center  of  the  mo.lcculc.  Our  geoinetry  differs  slightly  from  tjiat  used  in  ref.  13;  this. 
ct>iiki  be  the  cause  of  the  diflo'cnt  SO.MOs.  and  would  also  aflcct  ihcg-tensor  calculations  (sec  below). 

,  .  .  .  More  than  half  of  the  SD  in  Ga^As  resides  on  the  As  (0.5 1  5),  withO.  242  per  Ga  atom.  A  similar  trend, 
is  rep  oiled  in  ref  1  3,  but  with  more  SD  on  the  apical  .As  atom:  0.694  per  As  and  0.  153  perGa,  with,  larger 
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{inisotropic  than  isotropic  SD  contributions  for  each  atom  A;,,  is  largest  for  the  Ga  atoms,  and  the  spin-dipole 
couplings  are  smaller  for  Ga  than  As,  similar  to  the  results  from  ref.  13. 

4.3.  Ga^AS)  '  .  ’  .  ,  ■  ■  ■  '  '  ■  ■  •  • 

Our.  results  place  the  majority  of  the  SD  on  the  Ga  atoms'  (0.369),  with  only  0.088  per  As  atom.  This 
.  agrees  with  the  SD  distribution  from  ref  12  (0.314  per  Ga,  0.124.  per  As).  Based  on  the  ejqjerinental  results 
it  was  estimated  that  the  unpaired  electron  was  mainly  confined  to  p  a  orbitals  on  the  axial  galliumatoms.^  This 
estimation  is  not  supported  by  the  SOMO  (as  also  noted  in  ref  12),  which  shows  sizeable  in-phase  7t -bonding 
between  the  equatorial  As  atoms  via  p^  orbitals,  and  Po  antibonding  between  Ga  and  As  atoms. 

6ur’Ai,'„  differs  from  the  Cxperimental  "GajAs’j  (Ar  matrix)?  resiilfby  13%.  Ref.  12  repots  calculate 
'■  Aj,„’  varues  'for  only  ’''Ga,  even'thdugh  the  mtural  abundaricfe  is  60.4% ‘’Ga  and  3.9.6%  ’'Ga;  their  Aj,„' Value 
is  almost  identical  to  the  derived  experimental  ’’GujAs,  result  The  experimental  results  were  deriyed  assuming- 
.gii  =  &>  3nd  A||  s  A±  s  Aj,„  for.  As.-  This  does  not  appear  to  be  valid,  based  on  results  fi-om  ref  12  wliere 
anisotropic  spin  distributions  for  the  As  atoms  were  calculated  to  be  ca.  twice  the  isotropic  spin  distribution. 

4.4.  GajAs2  .  ,  .  \  •  •  .  ■  ' 

Tlie  small  distortion  from  D,,,  to  C,,  symmetry  makes  a  significant  impact -on  the  SD  and  charges  of 
theGa  atoms  in  G'ajAsj,  as  well  as  the  SOMO.  The  SOMO  is  conposedpf  in-phase  p^  orbitals  on  Ga,  and  the 
As'atom.s,  with  p„  orbitals  bn  Ga,  and  Ga/oriented  to  form  Ga.-As  and  Ga-As  antibonds.  The  SD  is  foiir 
times  larger,  and  the  charge  is  two  times  smalla-,  for  the  Ga,  atom  (Table  2)  than  for  Ga,  and  Ga,.  the' A,,,, 
.lor  Ga,  is  much  largci'  than  for  Ga,  and  Ga,,  and  tlte  SD,  charge,  .A,,.,,  and  anisotropic  (T,, hfcc  v.alues  for 
Ga,  arc  all  similar  lo  that  for  each  of  iheGa  atoms  in  Ga.As,. 


■i.5.  CiaASj 
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liie  SOMO  of  Ga As^  lies  in  the  yz-plane  from. orbitals  on  all  atoifis;  the  orbitals  on  Ga,  As,  and 
As^  are  in-phase,  and  of  opposite  phase  to  those  on  Asj  and  AS4. 

The  majority  of  SD  lies  on  As,  and  Asj  (see  Figure  1),  with  very  little  on  the  other  three  atoms.  This 
is  reflected  in  the  larger  A^i^  values  for  these  atoms  ^able  2),  whicfr  ^e  an  ordCT  of  magnW  larger  for  ,  - 
and  As,  than  for  As,  and  As.. 


4.6.  Ga^As 

The  SOMO  lies  in  the  plane  of  the  molecule,  with  p^  orbitals  of  Ga,  and  Ga^  oriented  to  be  bonding 
between  Ga.-Ga,,  with  the  po  orbitals  of  Ga4  and  Gaj  positiaied  to  be  Gai-Ga.,  and  Ga^-Ga^  antibonding.  The 

p;  orbital  on  As  is  small  arid  is  oriented  to  be  in  phase  with  the  py  orbitals  between  Ga,  and  Gaj. 

■  •  Them  is  alrriost  no  SD  associated  with  the  As  atom  iri  Ga4  As,  but  the  largest  a  tomic  charge  r  esides  on 

As  (-0.509).  Most  of  the  SD  resides  on  Ga,  and  Ga^,  and  is  dd.  twice  the  SD  for  Ga,  and  Ga4;  whereas  the 
charges  on  Ga,  and  Ga.  is  ca.  half  that  for  Gaj  and  Gaj.  This  invase  relatiaiship  between  charge  andSD  was 

obse^ed  by  Bruna  and  Greinforalkaii-metal  diatonic  radical’ cations.^"'  '  ■  .  •  • .. 


.  5.  G-tensors 

■  In  Table  3,  the  types  of  excited'  states  that  couple  with  four  different  ground  states  in-  C,„  symmetry 
(XA,,  X^A.,  XA,,  X’B,)  are  listed.  Table  4  summarizes  our  total  Ag  (2"‘‘  order)  values' and  compares  them 
to  known  experimental  and  theoretica’l  results.  To  oiir  knowledge  there’are  no  'experimental  or  calculated  ESR 
data  for  Ga,AS;,  GaAs^,  or  Ga'/A-i  For 'all  Ga.ASj.  doublet  radicals*(x  y  =  3,  5)  AE,  SO,  'L  and  Ag  values 
lor  the  first  five  excited  states  for  each  iireduciblc  are  given  in  Tables  5  -  10,  and  results  for  higher  excited 
siate.s  if  they  have  a  strong  maunelic  coupliiig  with  the  ground  stale.  Fifteen  excited  slates  were  calculated  for 
_  GujASj,  GaAs4  and  Ga^As  since  significant  magnetic  coupling  with  the  ground  slate  was  still  obseryed  in  tlie 
higher  staiesof  initial  9-root  calculations.  Tables  5-10  can  also  be  u.sed  fpr  vertical  excitation  ciiergies,  which 
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ai€given  inthebc%oftbLetableor  asa  foota^^  *  ' 

All  Cjv  molecules  have  the  z-axis  along  the  Cj  symmetry  axis.  GaAsj  and  Ga2As  are  placed  in  theyz- 
plane.  •  '  '  . 


5.7.  GaAs2  .  .  .  .  .  .  ■  • 

Table  5  shows  that  for.Ga  Asj  is  governed  by  the  .coupling, of  . li* A,  (.4a,  ->  SOMO-I  to 

SOMO)  and  2'A,  (3a,  -  2bj,  SOMO-3  to  SO  MO)  with  X'B^. 

The  Ag„,  component  is  dominated  by  the  coupling  with  1  -  2b,),  SOMO  to  LUMO+1 .  There 

are  also  strong  couplings  with  2"A,  and  3'A2  but  they  almost  cancel  each  other:  their  individual  Ag^^ 
cdntributioris  are  siMlar  in  magtiitude  and  opposite  in  sigh.  This  is  a  general  trend  dbservedfor  a  pair  of  states  ' 
gdierated  by  a  three  dperi  sh'ell-confi^ration; e.g.  for  HjCO^  the  A'E,  SO  arid  L  values  for  the  1,2^A2  ' 
•states  (Agyy)  were  of  similar  ihagnitude  but  their  contributions  to  Ag  were  of  opposite  sign.-^®  In  thecase  of 
GaAsj,  tlie  AE  and  SO  values  for  2- A,  and  3-A2  are  sirnilar,  but  the  L  values  differ  by  ca.  50%  (15  %  deviation 
in  L  values  in  the  HjCO*  example^’).  This  variation  of  the  L  values  may  occur  since  although  the  leading 
•configuration  for  both  23A2-and  3^A'-  is  a  4a,  -  2b,-  excitation  (90%  and.  8.5%,  respectively),  .there  are  small  . 
differences  iri  tlie  other  configurations  contributing  to  2^A2  and  3^A2  (Table  5). 

The  total  Ag„  component  is  an  order  of  magnitude  smaller  than  Ag;;  and  Ag„,.  Althaigh  thestrongest- 
couplings  withX’Bj  arise  from  5‘B,  and  6^B|,  theirAg  contributions  almost  cancel  each  other,  since  both  of 
these  states  derive  from  a  three  open  shell  configuration’'^^'’''''’  (3a,  -  la';.,  SOMO-3  to  LUMO).  In  this  case  the 
AE  for  tlie  two  states  are  very  close,  but  the  SO  and  L  values  differ  by  cci'.  50%.  These  differences  arise  from 
mixing  of  the  leading  configuration. 

Calculated  g-sliilts  foi  GaAs,  by  Arratia-Pcrcz  et  «/'■’  utilized  a  lirst-ordcr  perturbation  procedure 
based  on  the'SgF-DSW-X.a-  method,'^  and  a  geometry. by  .Lou  et  Altliough  tliiee  g-components  are 
expected  lor  C,,  symmeliy,  they  give  only  Ag,  and  Ag±  values.  Assumingtheir  Ag,  is  along  As-.As  (Ag^J  and 
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Agx  IS  Ag,„  our  calculated  Ag  values  for  GaAsj  differ  significantly  fixun  those  calculated  in  ref  13,  being 
ca.  70000  ppm  larger  in  magnitude  There  is  no  experimental  ESR  data  available  for  GaASj  to  gauge  the 
accuracy  of  dther  result.  ' 

S.lGa^As  •  •  •  •  .  .  ’ 

A.g,,;  value  of  GajAs  is  not  dominated.by  coupling  widi  any  one  excited  state;  rather  PA,  -  4  ^A, 
all  make  large  negative  contributions  to  Ag„,  countered  to  a  small  extent  by  6^A,  (Table  6). 

The  two  largest  contributions  to  Ag^y  are  negative  values  fiom  2^A2  and  4^A2,  countered  by  a  positive 
one  from  S’Aj.  As  die  leading  configurations  for  4^ A,  and  S'A^  are  the  same  (3a,  -*  2b|),  the  AE  and  L  values 
are  Similar  and  thdr  A  g  contributions  are  of  opposite  sigh, 'but  the  SO  values  differby  cn.'50%.  This  is  due 
to  S^Aj  being  dominated  by  a  three  Open- shell  configuration  (3  a,  ^  2b,), 'Whereas  4^A2  is 'more  evenly' 
'comprised  of  this  (40  %.)  and  a  double  excitation  from  the  SOMO-1  to  the  SOMO  and-LUMOt-l  (3a, ^ 
20^1  a,)  (37%).  ; 

The  Ag,2  value  is  dominated  by  coupling  to  PB,,  with  contributions  from  3,5,6^B,.TKe5^B,  andfr^B, 
states  have  the  Same  leading  three  open  shell  configuration,  SOMO*  1  (3a,)  to  theLUMO+2  (la^);  but  again 
the  A g. contributions  fr-om  tliese  two  states  do  not  cancel^®-’’-''"  due  to  diffa-ences  in  the  configuration  setup. 

As  with  GaAs,,  only  Ag.|  and  Agx  values  for  Ga,  As  are  given  in  ref  13.  If  we  assume  that  thdr  Ag,, 
is  along  Ga-Ga  ()'-axis,  ds'for  GaAsj  along  As-As),  then  the  Ag^y  conparison  is  vay'poor.  If  thdr  Ag,, 
'corresponds  iristead  to  Ag,„  then  the  result  differs  from  oiirs  by  ca.  20%,  but  is  incohsisteiit  with  the  GaAs", 
x,v,z-axis  assignments.  '  '  '  '  ’ 

5.3.  GaAs, 

In  Dj,;  syminetry  fX'Aj"),  l  lie  standardorieniaiion  places  the  G  a  atoms  ofGajASj  along  the  z-ax  is,  with 
the  As  atoins  in  the  .vr-piane;  .v  and  v  (and  Ag^,  and  Ag,.,)  arc  degenerate  (MZ')  Our  Ag  calculations  were  done 
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with  GajAsj  in  Cj,  symmetry  (X^B,),  placing  the  Ga  atoms  along  the  jc-tois  land  the  As  atoms  in  the^-plane 
with  one  As  atom  lying  on  the  z-axis.  In  this  case,  the  ’A,  and  states  correspond  to  ^E’.  For  checking  the 
accuracy  of  our  calculations,  we  calculated  Ag  for  both  ^A,  and  states.  The  degenerate  conponents  of  Ag 
will  always'  be  called  Ag,;  and  Agy^,  corresponding  to  the  notation.- Bue  to  the  indep^ident  seiectioil  of 
reference  configurations  and  ©ctrapdatipn  for  ^A,  and  ^Bj  states,  slightly  different  values  were  obtained  for 
Ag^andAg^j,.  . 

The  degree  of  degeneracy  of  Ag,,  and  Ag^^  can  be  checked  in  Table  7.  The  values  are  dominated  by 
coupling  of  the  ground  state  with  3‘A,  and  2^63,  respectively,  corresponding  to  excitation  from  the  SOMO 
(3b,)  to  the  degenerate  LUMO  (da,)  for  Ag„  and  LUMO+l  (3b2)  for  Agj.^.  The  second  largest  contribution 
is  from  5^A,,  an  exdtatibn  fran  the  SOMO  to  the  degenerate  LlMO+4  (7a,),  and  from  4^63,  SOMO  to 
LUMO+3  (Sbj).  The  remaining  maja  excited  S'tate  cbuiilings  iiivolve  b.f'Aj  and  5, d^B^,  excitations  from  the  '' 
degenerate  SOMO-2  (la2,'Ag,^)  and  SOMO-1  (Tb,,  Agzz)  to  the-LUMO+2  (4bi),  resulting  in  three  open  shell 
configurations.  As  shown  in  Table  7  for  Ag^^,  the  AE,  SO  a  nd  L  values  for  S'Bj  and  d’Bj  are  veiy  similar  and 
the  Ag  values  are  qjpositein  sign.  This  behavfour  is  ejpectedfoir  pairs  of  states  generated  by  a  three  open  shell 
configuration,^*-”'^“^This  is  not  as  evident  in.d^A,  and7^A|  in  Ag,„  where  theSO  and  L  values  are  similar,  • 
but  tlie  A  E.  values  differ  by  .  1.1  eV.  •  .  .  .  •  • 

As  shown  in  Table  7,  Ag„  was  calculated  to  be  an  order  of  magnitude  smaller  than  Ag^^  and  Ag„,' 
with  the  largest  contribution  from  coupling  with  the  3’ A^  state,  an  excitation  from  SOMO-5  (43,)  to  LUMO+2 
(4^).  ■  ■  •  ••  •  .•••.•.: 

Ga^As,  is  the  only  Ga,As,  doublet  radical  w'ith  published  experimental  ESR  data;  comparison' of  our 
rcsullsto  tlic  experimental"  and  publishal  thoarctical'’  g-tensordata  is  given  in  Table  4.  Our  results  arc  in  good 
agreement  with  experiment,-  beinga  much  clo.ser  fit  than  the  theoretical  results  from  ref.  12. 

5.-/. 
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In  ovtr  calculations  tbree  Ga  atoms.pf  GajAsj  were  placed  in  the >^-plane  and  the  two' As  atoms 
along  the  x-axis.  The  distortion  from  aDj^  trigonal  bipyramid  to  Cj,  symmetry  was  accomplished  .by  having 
the  Ga,-Ga2  (and  Gai-Gaj)  distance  shorter  than  Gaj-Gaj  [See  Figure  l.Cj^Ca)] .  Ga,  lies  al  ong  the  z-axis,Ga2 
•and  Ga,  along y.  '  ■  ■''r' •  ' 

The  Ag„  value  is  governed  by  the  c.eiq)ling  of  l^Bj'and  3^62  with  the  X^A,  ground  state  (Table  8). , 
There  are  numerous  contributipns  to  Ag^y,  the  largest  being  l^B,  (,2b|  -  Sa,,  SOM  0-3  to  SO  MO)  and 
3^B|  (Sa,  -*  3b,,  SOMO  to  LUMO+2),  which  almost  cancel  each  other.  The  overall  Ag^y  is  relatively  small, 
with  most  of  the  small  negative  contributions  (2,4,6, 14^B,)  being  negated  by  positive  contributions  from  5^B, 
and  IPB,. 

The  ovetall  Ag^,  is  also  small,  dufe  to  cancdling  contributfons  from  pairs  of  states  generated  by  three 
'open  shell  configurations.  Contributions  to  Agj^  from  2*A2  and  3^A2,  both  having  the  same  configuration', 
almost  cancel  each  otter?  thdr  AE  values  arealmost  identfcal,  and  although  the  SO  of  2^A2  is  ca.  70%  larger 
than  that  of  3"  A2,  the  L  value  for  2’ A,  is  ca.  80%  smaller  than  tliat  of  S’Aj.The  leading  configuration  for  4'A2 
and  5.^A2  is  thesam'e  (laj'^  6a,,  SOMO-2  toLUMOfl),  but  the  resptetive  Ag’s  differ  vastly  in  magnitude, 
due  probabfy  to  some  mixing  of  configurations.  i  f  . 

'  5.5.  GaAs/  • 

The  molecule  was  oriented  such  that  Ga  lies  on  the  z-axis,  and  As,-As2  in  the  j'z-plane;  ASj  and  AS4 
lie  in  the  xZ-pl'ane  [Figure  1,  C,,  (a)].  ■'  •  • 

For  Ag^,,  magnetic  coupling  was  dominated  by  coupling  with  l^A,  (36,  -  63,,  SOMO  to  LUMO). 
Large  ncgtUive  contributions  from  2,5‘A|  arc  effectively  cancelled  by  4,9^^,  (Table  9). 

In  contrast  to  Ag^^,  the  overall  Agy,  and  Ag^^  values  were  small  due  to  cancelling  positive  and  negative 
Ag  components  from  numerous  excited  states^  For  Agyy,  large  contributions  from  1  ’A,  and2‘A2  are  ofqiposite 
(three  open  shells)  and  nearly  cancel.  The  next  largest  contributions  arise  from  5'A,  and  10' .A,,  both  being 
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positive.  However,  tte  overall  A  ^  is  relatively  small  and  negative  due  to  negative  contributions  from  a  number 
of  states  (3 ,4, 6, 7,9,1 1'Aj).  '  ' 

The  largest  contributions  to  from  PB,  and  2^B„  are  again  of  opposite  sign  although  the  two 
states  do  not  d^ive  from  ^e  same'con figuration.  largest  in  ma^itude  ^  from  S’B,  and 

1  I^B,,  which  have  the  sairie  three  open  shell  configuration.  Due  to  other  contributing  states,  the  final  Ag,^  is 
again  relatively  small. 

5.6.  Ga^A.<i 

In  our  calculations  Ga4 As  was  in  the  yz-plane,  with  the  z-axis  bisecting  the  Ga^-As-Caj  and  Ga4-As- 

Gaj  angles  .(see  Figure  1);  the  x-axis  was' perpendicular  to  the  molecular  plane.  ■  '  ■  . 

For  Ag„,  the  magnetic’  •couplihg  was  dominated  by '  two '  excited  states  (1,2^'2);  the  leading 
configurations  for  which  Were  the-same  (5a,  ^  4b2  and  Sbi  5a])f:but  contributed  in  an  approximate^  opposite 
manner  (I’Bj:  62%  and29%;  2-B2:  29%  and  63%, respectively,  seeTable  10). 

The  Ingest  contributions  to  Ag^^  were  from' hi^i-lying  excited' states  (6,7,8'.B,).  All  of  the  states 
•  contributing  significantly  to  Ag-^.  had  negative  Ag  components  with  the  exception  of  S^B,,  which  is  one  of  a 
pair  of  states  (-with  6^6,)  from  a  tliree  open  shell  configuration  (4ai  2b,). 

High-lying  states  also  dominate  Ag,.^,  with ‘the  largest  negative  contributions  from  6,7, 1 2,14^A2.  The 
largest  positive  Ag  contribution  is  from  1 1^  A2,  which  has  the  same  tliree  open  shell  configurat  ion  as  1 2- A,  (db, 
-  3b|),  although  neither  state  is 'dominated' by  this  configuration  (42%  for  1  PAj,  60%  for  1 2^  A,).  The  other 
state  with  a  positive  contribution  to  Ag„  is  4’ A,,  which  has  the  same  leading  three  open’shd'l  'configuration  as 
6-A,  (311,  -  2b|);  here  the  SO  values  are  quite  different  .and  the  AE’s  differ  by  ca.  1  eV,  and  two  Ag 
contributions  are  opposite  in  sign  but  are  not  .similar  in  magnitude. 


6.  Summary  and  Conclusions 
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New  programs  developed  for  the  calculation  of  g-tensors,  based  on  naultireference  Cl  wavefiinctbns, 
have  been  capable  of  detaming  g-tensors  for  III-V  doublet  radicals  with  up  to  163  electrons  in  an  efficient 
manner,  and  without  the  need  for  RECPs  (although  this  may  still  be  necessary  for  the  study  of  larger  clusters). 

In  contras  t  to  smaller  molecules  for  which*  strohg  magnetic  coupling  with  the  ^und  state  is  observed 
with  on\y  a  few  excited  states  (e.g.  refs.  28-30),  signific^t  contributions  to  Ag  were  often  observed  .from 
various  higher  excited  states  of  Ga^ASy.  This  is  due  in  part  to  the  higher  enagy-density  of  excited  states  of 
these  molecules,  e.g.  14  excited  states  are  within4  eV  for  Ga4As  (Ag^^  in  Table  10) . 

From  comparison  with  experimental  ES  R  data  for  Ga.  AS3  ^  (Table4)  it  is  evident  that  our  new  g-tensor 
results  are  significantly  better  than  the  DFT  results  reported  in  ref.  12,  although  they  coixectly  reproduced  the 
‘  trends  of  the  experimental  results  (/;e  ’Ag^^  «A  gyy,  Ag^^).  This  suggests*that  our  g-tensof  results  for  GaAs2  ' 
and  GajAs'arfe  also  improved  Over  pubHshedDFT  values, and  that  results:  for  the  other  Clusters  are  reliable. 
Future  work  includes  using  these  new  programs  investigate  the  structures  and  properties  of  larger  III-V  • 
clusters  (1 1  + atoms),  and  tte  study  of  paramagnetic  defects  (heteroatomic/  anti-site^'^  and  hoiedefects‘0  within 
these  larger  clusters.  Results  presented  m  this  paper  show  that  one  does  not  need  to  consider  all  valence 
electrons  (VE),  which  decreases  calculation  times  and  may  allow  larger  molecules  to  be  calculated  before  the 
use  of  I^CPs  is  required.  .  .  •  . 

Also  planned  are  modifications  to  the  program  code  to  allow  g-tensor  calculations  for  triplet  and  ’ 
quartet  multiplicities.  AnumberofAB' III-V  compounds  for  which  experimental  ESRdata  areavailable  (e.g. 
GaAs'^,^  GaP"  have  quartet  ground  states.  Higher  spin  states  are  not  restricted  to  ions,  e.g.'GaPj  has  a  ^A. 
ground  state."**  These  would  be  helpful  in  gauging  the  quality  of  our  results,  since  to  our  knowledge  Ga2As3  ^ 
and  BNB*  are  the  only  III-V  doublet  radicals  with  published  ESR  dita. 
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Table  1.  Optiriiized  bond  distances  (A)  and  an^es  (deg)  from  this  work  [KLYP/6-3H+<j(2dQ],  and 
comparison  to  literature  values,  for  all  Ga^As,  (x  +  y  =  3,  5)  doublet  radicals. 


Molecule,  symmetry,  ground  state 

Results 

GaAsj  Cj,, 

This  work 

■  Ref.  17  (MRCI) 

Ref  14  (MRCI) 

Ga-As 

.  2.775 

2.80  . .  • 

2.86 

As-Ga-As 

46.5 

45.9 

46.6 

GajAs  Cj,,  X% 

Thiswoik  . 

Ref.  17  (MRCI) 

.  Ref  19  (LSD) 

Ga-As 

2.384 

2.407 

2.33 

Ga-As-Ga 

96.3 

79.9 

98.0 

Ga^As,  . 

.  Thiswoik 

Ref.  25  (MRCI) 

Ref  19  (LSD) 

Ga-Ga 

4.268  •  • 

not  reported 

not  reported 

Ga~As 

2.594 

2.589 

2.65 

As-As 

2.555 

2.563 

•  2.62 

Ga3As2  C2v)'X^Aj 

Thiswoik  ' 

Ref  25  (MRCI) 

Ref  19  (LSD) 

Ga,“Ga2,  Gai-Ga3" 

.3.751 

3.702 

3.72 

Ga2-Ga3 

3.979 

4.114  ■ 

4.57 

Ga,-Asi 

2.446 

2.401 

2.41 

Ga2-Asj,  Ga3“As, 

2.67 1  ■ 

.  2.725 

2.59. 

Asi-As2 

2.725 

2.782 

2.70 

Ga2-Ga’|-Ga3 

64.0 

•67.5 

70.0 

AS|-Gai-As2 

67.7 

70.8 

67.9 

ASi-Ga2“As2,  As,-Ga3-As2 

61.3 

61.4 

55.2 

Table  1  continued. 


GaAs4  Cjv,  X^B2 
GaTASj,  Ga-Asj" 

Ga-Asj,  Ga-As4 

Asj-Asj,  AS1-AS4,  AS2-AS3,  As2-As, 
AS1-AS2 

Ga4As  Cjv,  X^A, 

Ga,-As,  Ga2~As" 

Gaj-As,  Ga4-As 
G^aj“Cfa3j  G^a2“Ga4 
Gai-Ga2 

Ga,-Ga4,  Ga2-Ga3 

Ga3-Ga4 

Gaj-As-Ga2 

Ga,-As-Ga35  Ga2-As-Ga4 
Gai-As-Ga4,  Ga2-As-Ga3 
Ga^-As-Ga. 


This  woric 

Ref.  15 

3:116 

no  data  given: 

2.545 

C2V,  edge-capped 

2.454 

tetrahedron 

3.006 

This  woik 

Ref  15 

2.588 

no  data  given: 

2.487 

C2V,  planar 

2.870 

.  trapezoid 

2.599 

with  As  in  center 

4.583  ‘ 
4.913 
60.3 
68'.8  ; 
129.1 
162.0 


a  See  Figure  1  for  atom  labelling  scheme. 
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Table  2.  Spin  densities,  atomic  charges,  and  hyperfine  coupling  constants  (MHz)  for  all  Ga^As,  (x  y  =  3, 5) 
doublet  radicals. 


Molecule 

Atom 

Spin  Density 

Charge" 

A.„ 

T„/T.„/T./ 

GaASj 

Ga 

0:262 

0.226  . 

26 

-57/ 131/ -74 

Ref  10 

0.338 

0.193 

83“' 

As 

0.368 

-0.M3 

7 

••  -lll/-113'/225' 

Ref  10 

0.331 

-0.096 

125“' 

GajAs  Cjv, 

Ga 

0.242 

0.128 

447 

-63/ -81'/ -63' 

Ref  10 

0.153 

0.103 

ir- 

As 

0.515 

-0.257 

-148 

-159/ 304/ -145 

Ref  10 

0.694 

-0.206 

216'' 

Ga^As,  D3„  X^Aj" 

.  Ga 

■  0.369 

0.194 

1325" 

-65  /  -65  /  129 

Ref  9  . 

.0.314 

0.150 

.  sr" 

As 

0.088 

-0.129  • 

•  -23" 

-35/-36/71 

Ref  9 

0.124 

-0.100 

14" 

GajAs, 

Ga^ 

0.426 

0.175 

1735 

-85/ -71  /  156 

Gaj,  Ga3 

0.085 

0.276 

203 

■■  -21/ 45'/ -24' 

ASi,  ASj 

0.202  • 

•  -0.363 

-25 

-61'/ -56/ 116' 

GaAs^  C2v,X'B3 

Ga 

-0.038 

0.214 

-lb 

4/ 17/ -21 

Asj,  Asj 

0.585 

-0.078 

-46 

-156/ 308'/ -152' 

ASjj  A-S^ 

-0.066 

-o;o3o 

12 

-10'/ 45/ -36' 

Ga,As  Q^X'A, 

Gaj,  Ga2 

0.342- 

0.083 

-122 

•  -89/  177'/ -88' 

Ga3,  Ga^ 

0.169 

0.171 

153 

-34  /  -43'  /  76' 

As 

-0.021 

-0.509 

0.05 

.  -13/ -7/ 20 

a  Atomic  charges  from  a  Mulliken  populatiai  analysis. 

b  Off-diagonal  spin-dipolar  (anisotropic)  coupling  terms  T'„,  T',,  are  indicated  by  a '  after  tlie  value. 

(■'  Ai,„=  1524  [1936]  MHz,  Aj^^  =  87  [71]  MHz  from  experimental  ESR  spectra  for  ^’GajAsj  f  Ga;,As3].^  A^,^, 
=  2012  MItz,  Aj;p  =  88  Mliz  for  ’’Gaj^Asj  in  ref  12. 
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rflbe  spin-dpolar  contribution  to hfccfrdmref.  13,  defined  as  the  expectation  value  (3x^-1^)/!^  averaged  over 
the  relativistic  wave  function. 
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Table  3 .  Excited  state  couplings  for  the  four  possible  ground  states  in  Q,  symmetiy. 


Ground  State 

Ag„ 

Excited  State  Coupling 

Agyy 

> 

OQ 

N 

N 

X^A, 

'B,  ’ 

/A, 

X'Aj 

X'B, 

^A, 

X'Bj 

^A, 

’B, 
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Table  4.  Calculated  g-tensor  data  (Ag  in  ppm)  for  Ga^As^  (x  +  y  =  3,  5),  and  comparison  with  experimental 
and  other  theoretical  results. 


Agxx 

Ag„ 

A&. 

GaAsjCX^BO 

175300 

-175120 

-18280 

Ref.  13,  calc. 

95600 

-116500 

Ga^As  (X'Bj) 

-125950 

-24850 

51030 

Ref  13,  calc. 

-188500 

40100 

GajAs3(X^A2"inD3j 

-73410“ 

- 

6460 

Ref.  2,  expt. 

-82300 

- 

~0 

Ref  12,  calc. 

-148000 

- 

-16200 

Ga3Asj  (X'A,) 

-71590  ■ 

-11150 

■  13270 

GaAs4  (X'Bj) 

-171030 

-7540 

-14720 

Ga^As  (X'A,) 

-123540 

-21850 

-43470 

a  Average  of  Ag,,  and  Ag^y  from  Table  7. 


Table  5.  Calculated  values  of  AE,  SO,  L  and  Ag  (2"'*  order)  for  excited  states  having  a  large  magnetic  coupling 
with  the  ground  state  of  GaAsj  and  vertical  excitation  energies”. 


Component 

AE(eV) 

SO  (cm‘) 

L(au) 

Ag  (ppm)* 

Ag„eA,) 

PA,  (4a,  -  2b2,  S-1  -  S)” 

1.47 

169.5 

1.096 

64366 

2'A,  (3a,  -  2b2,  S-3  -  S) 

2.67 

475.1 

1.105 

99952 

3'A,  (lb,  -  la2,S-2  -  L) 

2.89 

165.1 

0.149 

4333 

4^A,  (2b2  -  5a„  S  -  L+3) 

3.97 

116.2 

0.784 

11687 

5' A,  (lb,  -  la^,  S-2  -*  L)  [30%] 

4.24 

69.2 

-0.404 

-3357  175300 

Agyy  CAj) 

PAj  (2b2  -  laj,  S  -  L) 

1.26 

476.4 

-0.897 

-172824 

2'Aj(4a,  -  2b„S-l  -.L+1) 

3.19 

98.4 

0.872 

13686 

3'A2  (4a,  -  2b„  S-1  -  L+1) 

3.90 

-86.1 

1.530 

-17199 

4^Aj  (3a,4a,  -  2bjlaj,  S-3  +  S-1  -  S  +  L) 

4.37 

31.8 

0.034 

128 

5'A2  (4a,"  -  2b2la2,  S-1  -*  S  +  L)  [26%] 

4.71 

50.1 

0.130 

708  -175120 

AgzzCB.) 

PB,  (2bj  -  2b„  S  -  L+1) 

2.27 

-129.3 

0.657 

-19067 

2"B,  (lb,  -  2b2,  S-2  -  S) 

2.55 

36.5 

-0.042 

-307 

3"B,  (4a,  -  laj,  S-1  -  L) 

2.62 

81.2 

-0.033 

-526 

4"B,  (4a,  -  la^,  S-1  -  L)  [57%] 

3.13 

250.8 

-0.149 

-6072 

5"B,  (3a,  -*  laj,  S-3  L) 

3.61 

348.0 

0.983 

48315 

6"B,  (3a,  -  la^,  S-3  -  L) 

3.89 

190.9 

-1.612 

-40264  -18280 

a  The  ground  state  is  ...lb2Mb|Ma|^2b2'  (13  VE).  Vertical  excitatiai  energies  for2^Bj  to  S^Bj  are  3.58, 
4.25,  4.35,  and  4.76  eV,  respectively. 


b  Total  contribution  for  all  calculated  excited  states  in  boldface. 


cS  =  SOMO,L  =  LUMO. 
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Table  6.Calculated  values  of  AE,  SO,  Land  Ag  (2"'’  order)  for  excited  states  having  a  large  magnetic  coupling 
with  the  ground  state  of  GajAs  (X^Bj),  and  vertical  excitation  energies". 


Component 

AE(eV) 

SO  (cm  ') 

L(au) 

Ag  (ppm)‘ 

Ag„  CA.) 

PA,  (3a,  -  2bj,  S-1  -  S)' 

0.17 

177.4 

-0.069 

-36004 

2'A,  (2b2  -*  4a„  S  -*  L)  [58%] 

1.97 

243.5 

-1.019 

-64115 

3^A,  (lb,  -  la^,  S-2  -*  L+1)  [33%] 

2.72 

-60.7 

0.735 

-8368 

4^ A,  (3a,'  -  2bj4a„  S-1  -*  S  +  L) 

2.90 

-204.4 

0.910 

-32678 

5'A,  (lb,  -  laj,  S-2  -  L+1) 

3.38 

-30.86 

0.306 

-1424 

6'A,  (3a,  -  3b2,  S-1  -  L+3)  [39%] 

3.72 

84.9 

1.271 

14761 

-125950 

Agyy  CAy) 

I'Aj  (lb,  -  4a„  S-2  -  L) 

2.32 

2.7 

0.154 

91 

2'A2  (2b2  -  laj,  S  -  L+1) 

2.57 

156.8 

-0.477 

-14834 

3'A2  (lb,  -  4a„  S-2  -  L) 

2.78 

41.9 

0.014 

108 

4'Aj  (3a,  -  2b„  S-1  -  L+2)  [39%] 

3.32 

158.5 

-0.784 

-19071 

S'Aj  (3a,  -  2b„  S-1  -  L+2) 

3.82 

100.8 

0.737 

9909 

-24850 

Ag..eB,) 

PB,  (lb,  -  2b2,  S-2  -  S) 

0.33 

254.7 

0.219 

87664 

2'B,  (lb, 3a,  -  2b24a„  S-2  +  S-1  -  S  +  L) 

2.57 

-0.2 

0.140 

-4 

3'B,  (3a,  -  laj,  S-1  -  L+2) 

2.66 

122.2 

0.904 

21148 

4'B,  (lb,3a,  -  2b24a„  S-2  +  S-1  -  S  +  L) 

3.01 

41.3 

0.051 

360 

5'B,  (3a,  -  laj,  S-1  -  L+2) 

3.32 

-220.1 

1.066 

-35932 

6'B,  (2b2  -  2b„  S  -  L+2)  [40%] 

3.74 

146.9 

-1.052 

-21056 

51030 

a  The  X^Bj  ground  state  is  ...lb2^1b,^3a,^2bj'  (11  VE).  Vertical  excitaticxi  aiergies  for2^B2  to  S^Bj  are  2.24, 


2.65,  3.17,  and  3.41  eV,  respectively. 
b  Total  contribution  for  all  calculated  excited  states  in  boldface. 


cS  =  SOMO,L  =  LUMO. 
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Table  7.  Calculated  values  ofAE,  SO,  Land  Ag(2"‘‘ order)  for  excited  states  having  a  large  magnetic  coupling 
with  the  ground  state  ofGajAsj  and  vertical  excitation  energies". 


Component 

AE(eV) 

SO  (cm  ') 

L(au) 

Ag„  CA,) 

PA,  (5a,  -  3b„  S-3  -  S)" 

1.98 

41.5 

-0.037 

2'A,  (4a,  -  3b„  S-5  -  S) 

2.98 

0.4 

-0.014 

3'A,  (3b,  -  6a„  S  -*  L) 

3.30 

342.9 

-1.777 

4' A,  (3a,  -  3b„  S-6  -*  S) 

3.87 

0.02 

0.003 

5'A,  (3b,  -  7a„  S  -  L+4) 

3.70 

246.7 

0.634 

6'A,  (laj  -  4b2,  S-2  -  L+2) 

3.43 

148.1 

0.396 

7' A,  (laj  -  4bj,  S-2  -L+2) 

2.29 

133.4 

-0.368 

Ag,,(^B,) 

PBj  (2bj  -  3b„  S-4  -  S) 

2.69 

-40.5 

0.037 

2^Bj  (3b,  -  3b2,  S  -  L+1) 

3.28 

-355.5 

1.761 

3^82  (3b,  -  4b2,  S  -  L+2) 

3.50 

4.9 

0.032 

4^Bj  (3b,  -  5bj,  S  -  L+3) 

3.67 

269.0 

0.713 

5^82  (2b,  -  4b2,  S-1  -  L+2) 

3.83 

101.8 

0.316 

b'Bj  (2b,  -  4b2,  S-1  -  L+2) 

3.93 

119.2 

-0.342 

Ag  (ppm)* 


-398 

-1 

-94074 

0 

21549 

8742 

-10911  -75230 

-284 

-97100 

23 

26635 

4287 


-5280  -71590 
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Ag„eA,) 


l^Aj  (laj  -  3b,.  S-2  -  S) 

1.58 

-0.1 

3x10^ 

0 

2*A2  (5a,  -  4b2,  S-3  -  L+2) 

3.64 

0.9 

6x10’ 

0 

3'A2  (4a,  -  4b2,  S-5  -  L+2) 

2.89 

71.8 

0.401 

5083 

4^A2  (3b,  -  2a2,  S  -  L+7) 

4.55 

0.3 

2x10^ 

0 

5’A2  (5a,  -  4b2,  S-3  -  L+2) 

3.85 

0.7 

-0.021 

-2 

b'Aj  (5a,  -  5b2,  S-3  -*  L+3)  [50%) 

4.09 

194.5 

0.084 

2041 

a  The  grotind state  is  ...2b2^5a|^la2^2b|^3b,'  (21  VE).  Vertical  exdtation  energies  for2^B,  to  5^B,  are 
1.58,  3.85,  4.54,  and  4.50  eV,  respectively. 
b  Total  contribution  for  all  calculated  excited  states  in  boldface. 


c  S  =  SOMO,  L  =  LUMO. 


Table  8.  Calculated  values  of  AE,  SO,  L  and  Ag  (2"'’  order)  for  excited  states  having  a  large  magnetic  coupling 
with  the  ground  state  ofGajAsj  (X^A,),  and  *A,  vertical  excitation  energies". 


Component 

AE(eV) 

SO  (cm-') 

L(au) 

Ag  (ppm)* 

Ag„CB,) 

VB,  (2bj  -  5a„  S-1  -  S)" 

0.72 

-72.2 

0.525 

-26906 

2^Bj  (la^  -  3b„  S-2  -  L+2) 

1.87 

5.1 

0.038 

53 

3^6^  (5a,  -*  3bj,  S  -  L)  [48%] 

2.94 

221.4 

-1.402 

-53718 

4'Bj  (laj  -  3b„  S-2  -  L+2) 

2.52 

46.7 

-0.653 

-6177 

5'Bj  (2bj.-  6a„  S-1  -  L+1)  [52%] 

2.94 

22.2 

-0.921 

-3537 

6'Bj  (2b2  -  6a„  S-1  -  L+1) 

3.16 

72.1 

1.605 

18664  -71590 

Ag,,(^B,) 

PB,  (2b,  -  5a„  S-3  S) 

1.15 

80.3 

0.902 

32076 

2'B,  (laj  -  3bj,  S-2  -  L) 

2.52 

41.2 

-0.783 

-6518 

3'B,  (5a,  -  3b„  S  -  L+2) 

2.94 

320.5 

-0.599 

-33328 

4'B,  (laj  -  3b2,  S-2  -  L) 

2.58 

15.3 

-0.506 

-1530 

5'B,  (2b,  -  6a„  S-3  -  L+1) 

3.90 

24.8 

0.567 

1840 

6'B,  (4a,  -  3b„  S-4  -  L+2) 

2.50 

-79.7 

0.511 

-8297 

ll'B,  (5a,  -  4b„  S  -  L+3) 

4.34 

71.2 

0.955 

7981 

14'B,  (3a,  -  3b„  S-5  L+2)  [26%] 

4.20 

73.5 

-0.732 

-6536  -11150 
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Ag,,  CA2) 

(1%  -  5a„  S-2  -  S)  1.00 

2'Aj  (2b2  -  3b„  S-1  -  L+2)  2.73 

3^A2  (2b2  -  3b„  S-1  -  L+2)  2.79 

4*Aj  (la^  -  6a„  S-2  -  L+1)  3.05 

S^A;  (la^  -  6a,.  S-2  -  L+1) _ ^ 

a  The  X^A,  ground  state  is  ...4a'^2b,^la2^2b2^5a,'  (19  VE).  Vertical  excitation  energies  for  2^A,  to  S’Aj  are 
1.68,  2.18,  2.56,  and  2.90  eV,  respectively. 
b  Total  contribution  for  all  calculated  excited  states  in  boldface. 
cS  =  SOMO,L  =  LUMO. 


76.5  0.394  15339 

253.8  0.252  11927 

75.3  -1.121  -15409 

20.8  0.630  2191 

2.7  -0.415  -200  13270 
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Table  9.  Calculated  values  of  AE,  SO,  L  and  Ag  (2"'’  order)  for  excited  states  having  a  large  magnetic  coupling 
with  the  groimd  state  of  GaAs4  and  ^82  vertical  excitation  energies. 

Component  AE  SO  L(au)  Ag(ppm)* 

_ _  (eV)  (cm‘) 


Ag„  CA.) 

I'A,  (Sbj  -  6a„  S  -  L)' 

2'A,  (5a,  -  3b2,  S-1  -  S) 

3'A,  (3b2  -  7a„  S  -  L+4) 

4'A,  (2b2  -  6a„  S-3  -  L) 

5'A,  (2b2  -  6a„  S-3  -  L)  [64%] 

9'A,  (5a,  -  4b2,  S  -  L+1)  [39%] 

Agyy  CAi) 

VA^  (3b,  -  6a„  S-2  -  L) 

2'A2  (3b,  -  6a„  S-2  -  L) 

3'A2  (3b2  -  2a2,  S  -  L+5) 

4'A2  (2b,  -  6a„  S-4  -  L) 

5'A2  (2b,  -  6a„  S-4  -  L) 
b'Aj  (5a,  -  4b„  S-1  -  L+3) 

T'Aj  (3b,  -  7a„  S-2  -  L+4) 

9'A2  (laj  -  3b2,  S-5  -  S) 

10^ Aj  (3b,  -  7a„  S-2  -  L+4) 
ll^Aj  (5a,  -  4b|,  S-1  -  L+3)  [55%] 

IS'Aj  (3b,3b2  -  6a,5b2,  S-2  +  S  -  L  +  L+2)  [28%] 


0.55 

390.3 

-0.496 

-179018 

2.22 

-249.0 

1.0471 

-59761 

2.46 

171.5 

0.070 

2501 

2.82 

265.0 

0.7923 

37977 

2.61 

-139.8 

0.386 

-10558 

4.10 

162.8 

1.850 

37462  -171030 

1.91 

-311.7 

0.807 

-66989 

2.32 

162.1 

1.378 

49104 

3.06 

-188.6 

0.117 

-3699 

3.37 

336.6 

0.057 

-2937 

3.18 

129.2 

0.794 

16423 

2.91 

123.4 

-0.372 

-8027 

3.44 

111.3 

-0.321 

-5292 

3.78 

49.5 

-0.355 

-2363 

4.43 

83.0 

1.561 

14898 

3.67 

84.2 

-0.494 

-5790 

3.66 

70.5 

0.702 

6890  -7540 

36 


Ag^CB.) 


I'B,  (3b,  -  3b2,  S-2  -  S) 

1.90 

285.7 

0.777 

59501 

2^B,  (3b2  -*  4b„  S  -  L+3) 

2.57 

-259.9 

1.285 

-66074 

3*B,  (3b,3b2  -  6a, ^  S-2  +  S  -  L) 

2.40 

163.7 

0.200 

6966 

4'B,  (2b,  -  3b2,  S-4  L) 

3.10 

293.5 

0.102 

4943 

5'B,  (3b,  -  4b2,  S-2  -  L+1) 

3.19 

-63.2 

0.765 

-7720 

6’B,  (3b2  -  5b„  S  -  L+6)  [58%] 

2.54 

-79.2 

0.921 

-14645 

8'B,  (3b,  -  5b2,  S-2  -  L+2) 

3.67 

-118.9 

1.378 

-22751 

ll^B,  (3b,  -  5b2,  S-2  -  L+2) 

3.72 

114.4 

1.319 

20686  -14720 

a  The  ground  state  is  ...Ia2^2b2^3b,^5a,^3b2'  (23  VE).  Vertical  excitatiai  energies  for  2*82  to  5^B2  are 
1.62,  2.60,  2.54,  and  2.63  eV,  respectivdy. 
b  Total  contribution  for  all  calculated  excited  states  in  boldface. 


cS  =  SOMO,L  =  LUMO. 


Table  10.  Calculated  values  of  AE,  SO,  L  and  Ag  (2""'  order)  for  excited  states  having  a  large  magnetic 
coupling  with  the  ground  state  ofGa4As  (X^A,),  and  *A,  vertical  excitation  energies". 


Component 

AE 

(eV) 

SO 

(cm-') 

L  (au) 

Ag  (ppm)‘ 

Ag„  CBi) 

l^Bj  (5a,  -  4hi,  S  -  L+1)"  [62%] 

1.28 

223.4 

-1.224 

-108538 

2*B2  (5a,  -  4bj,  S  -*  L+1)  [29%] 

1.59 

-51.6 

1.4029 

-23187 

(2b2  -  5a„  S-4  -  S) 

2.73 

-24.1 

1.423 

-6406 

4'B2  (lb,  -  lai,  S-3  -  L) 

2.98 

6.6 

-0.664 

-751 

5'Bj  (lb,  -  laj,  S-3  -  L)  [46%] 

3.02 

-23.0 

-0.147 

572  -123540 

Ag,yCB,) 

PB,  (5a,  -  2b„  S  -*  L+2) 

1.45 

-18.1 

0.216 

-1378 

2'B,  (lb,  -  5a„  S-3  -  S) 

2.23 

47.8 

-0.041 

-457 

3'B,  (3b2  -  laj,  S-1  -  L) 

2.03 

-1.6 

0.148 

-61 

4*B,  (3b2  laj,  S-1  -  L) 

2.44 

-30.0 

0.224 

-1405 

5'B,  (4a,  -  2b„  S-2  -  L+2) 

2.75 

249.0 

0.0782 

3601 

6'B,  (4a,  -  2b„  S-2  -  L+2) 

2.92 

-30.7 

1.616 

-8668 

7^B,  (5a,  -  3b„  S  -  L+3) 

3.06 

-72.1 

0.724 

-8675 

8'B,  (2b2  laj,  S-4  -  L) 

2.99 

-24.6 

1.409 

-5893  -21850 

T 

*  \> 


Ag^CA,) 

VA,  (5a,  -  laj,  S  -  L) 

1.22 

-63.9 

0.256 

31 

-6838 

2'A2  (4a,  -  laj,  S-2  -  L) 

2.49 

-0.06 

0.357 

-5 

3*A2  (4a,  -*  laj,  S-2  -*  L) 

2.62 

12.7 

-0.105 

-262 

4'Aj  (3b2  -  2b„  S-1  -  L+2) 

2.42 

-43.5 

-0.432 

3948 

S'Aj  (lb,  -  4b2,  S-3  -  L+1) 

2.93 

7.1 

-0.049 

-62 

6'Aj  (3b2  -  2b„  S-1  -  L+2) 

3.44 

215.6 

-0.422 

-13477 

7'A2  (5a,  -*  2a2,  S  -*  L+5)  [41%] 

3.26 

205.7 

-0.451 

-14511 

iPAj  (3b2  -  3b„  S-1  L+3)  [42%] 

3.83 

102.7 

0.847 

11592 

12'A2  (3b2  3b„  S-1  -  L+3)  [60%] 

3.96 

-63.9 

1.277 

-10487 

14* Aj  (2b2  -  2b„  S-4  -  L+2) 

3.98 

-39.4 

1.652 

-8334  -43470 

a  The  X*A,  ground  state  is  ...lb/4a,^3b2^5a,'  (17  VE).  Vertical  excitation  aiergies  for2^A,  to  5^A,  are  1.66, 
2.79,  2.80,  and  3.02  eV,  respectively. 

b  Total  contribution  for  all  calculated  excited  states  in  boldface. 


c  S  =  SOMO,  L  =  LUMO. 


Captions  for  ^nres: 


Figure  1 .  Optimized  geometries  of  Ga^ASy  (x  +  y  =  3, 5)  isomers,  shown  in  increasing  relative  energy  fromleft 
to  right.  Bond  lengths  and  angles  given  in  Table  1. 
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